arXiv:1506.08456vl [math.AP] 28Jun2015 


LINDA MARIA DE CAVEQ MARTA STRAN:0 


ASYMPTOTIC BEHAVIOR OF INTERFACE SOLUTIONS TO 
QUASILINEAR PARABOLIC EQUATIONS WITH NONLINEAR 

FORCING TERMS 


Abstract. We investigate the asymptotic behavior of solutions for quasilinear parabolic 
equations in bounded intervals. In particular, we are concerned with a special class of 
solutions, called interface solutions, which exhibit e metastable behavior, meaning that 
their convergence towards the asymptotic configuration of the system is exponentially 
slow. The key of our analysis is a linearization around an approximation of the steady 
state of the problem, and the reduction of the dynamics to a one-dimensional motion, 
describing the slow convergence of the interfaces towards the equilibrium. 


1. Introduction 


In this paper we study the asymptotic behavior of interface solutions to the initial¬ 
boundary value problem for quasilinear parabolic equations, i.e. 

dtu = edx {a{x)dxu) — G{u,dxu), x £ I = {—£,£), t>0, 

(1.1) au{±£,t) + bdxu{±£,t) = u±, t>0, 


_ u(x,0) = Mo(a:), 


X € L 


for some e,£ > 0, u± G M and a,b > 0. Concerning the function a{x), we require 
a G IE^’°°(I); as a consequence, there exist positive constants a, /3 G such that 

(1.2) a < a{x) < /3, for all x G / , 


so that the classical ellipticity and growth conditions are satisfied. Finally, concerning the 
initial datum uo and the nonlinear forcing term G, we assume 

uo G L^{I), G = Giz, w) G C^(m2). 

In particular, we focus our attention to the phenomenon known as metastability, whereby 
the time dependent solution develops into a layered function in a relatively short time 
(usually of order one), and then converges towards its stable configuration in a time scale 
that can be extremely long, depending on the size of the viscosity parameter e. 

Such behavior has been extensively studied for a large class of one-dimensional evolutive 
PDFs; to name some of these results, we recall here the area of viscous scalar conservation 
laws, with the contributions [llEKiailslIIIllSa], or phase transition problems, described 
by the Allen-Cahn and the Cahn-Hilliard equations equation (inEiEi ESI [miss]). 

Results on metastability for systems of scalar equations are less common; the slow 
motion for systems of conservation laws has been examined in m, while in [a [20] the 
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authors describe the phenomenon of metastability for systems with a gradient structure, 
with an analysis that is entirely based on energy methods. Finally, we quote [22], where 
the one dimensional Jin-Xin systems is analyzed. 

The bibliography is so rich that it would be impossible to mention everyone. 

Roughly speaking, the phenomenon of metastability can be summarized as follows: 
starting from an initial datum uq{x) which contains N zeroes inside the interval I, a 
layered solution with exactly N layers is formed in an 0(1) time scale; once this pattern 
is formed, it starts to move towards its asymptotic stable configuration, but this motion 
can be extremely slow as the viscosity parameter e goes to zero. As a consequence, we 
can distinguish two different time phases in the dynamics; a first transient phase of order 
one where the internal shock layers are formed, and a subsequent long time phase where 
the layers interact until the solution stabilizes to the stable steady state of the system. 

Usually, such behavior is related to the presence of a first small (with respect to e) 
eigenvalue of the linearized operator around the steady state (see, for example, [TO]). As 
it is well known, if Af is negative, the steady is stable; if, in addition, Af —?• 0 as e —)• 0, 
the steady state is metastable in the sense that the time dependent solution converges 
towards it in a time scale that goes to oo as e goes to zero. On the contrary, if Af is 
positive, the steady state in unstable and we will see the solution to ’’run away” towards 
a stable configuration; again, this motion will be extremely slow as e —)■ 0. 


The aim of this paper is to describe the metastable behavior of solutions to the general 
class of quasilinear parabolic equation described in (1.1). 

In the limit e —>• 0, equation (1.1) reduces to the hrst order hyperbolic equation 


(1.3) 


dtu = -G{u,dxu), 


complemented with initial datum uq{x) and appropriate boundary conditions. As it is 
well known, the set of solutions to (1.3) is the one given by the entropy formulation, in 
the sense of Kruzkov (see 


); moreover, the boundary conditions has to be interpreted 
in a nonclassical way in the sense of [3]. 

In the case a{x) = 1, for some special choices of the nonlinear forcing term G, it 
is possible to prove the existence of discontinuous stationary solutions for the inviscid 


problem (1.3), corresponding to stationary solutions with internal layers for the associated 


viscous problem (see, for example, m in the case of a reaction-convection equation); as 
already stressed before, the corresponding time dependent solutions exhibit a metastable 
behavior. 

There are a large number of works that have investigated such phenomenon for problem 
(1.1) with a{x) = 1; for instance, in [TO] [23[ |2l] j the authors describe this behavior for 
different parabolic equations, throughout the slow motion of the internal interfaces of the 
solutions. 


Motivated by this, we expect that also in the more general setting a{x) satisfying (1.2), 


all the discontinuities that appear at the hyperbolic level e = 0 shall eventually turn 
out into smooth internal layers, and that a metastable behavior will be observable in the 
vanishing viscosity limit. 
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In order to characterize the slow dynamics of solutions to (1.1), we mean to adapt the 


theory developed in m for parabolic equations to our more general setting; this strategy 
dates back the work of J. Carr and R. L. Pego [6] and can be summarized as follows. 

The principal idea is to construct a family i = 1,. ,N, of so called 

approximate steady states for the problem, and to linearize the original equation around 


an element of this family. With approximate steady state for (1.1), we refer to a solution 


that solves the stationary equation up to an error that is small in a sense that will be 
specified later. The parameters represent the location of the interfaces. 

The aim of this construction is to separate the two distinct phases of the dynamics. 
Firstly, we mean to understand what happens far from the steady state solution, when the 
interfaces are formed; subsequently, we want to follow the evolution of the layered solution 
towards the asymptotic limit. In particular, we describe such slow motion by obtaining an 

equation for the positions of the interfaces ^i{t), i = 1,. ,N. For the sake of simplicity, 

in this paper we restrict our analysis to the case of a solution with a single shock layer 
located in ^ £ I, the general case being similar (see, for example, 0). 

After the family {U^} is given, a study of the eigenvalue problem associated to the 
linearized operator obtained from the linearization is needed. Indeed, we want to 
describe the dynamics of solutions located far from the equilibrium conhguration of the 
system and, by performing a spectral analysis of C^, we are able to show that the speed 
rate of convergence of the solutions towards the asymptotic conhguration is small in e. 

We close this Introduction with an overview of the paper. In Section 2 we present 
the general strategy we develope in order to describe the long time behavior of solutions 
belonging to a neighborhood of a family of approximate steady states {U^{x]^)}^^j. By 


linearizing the original equation (1.1) around an element of the family, i.e. by looking for 
a solution u on the form u = + v, and by using an adapted version of the projection 

method in order to remove the growing components of the perturbation v, we obtain a cou¬ 
pled system for the variables (^, v), whose analysis is performed in the subsequent section. 
In particular, in Section 3, we state and prove Theorems 3.2[ 3.3 3.4 and |3.7[ providing 
different estimates for the perturbation v, depending on the choice of the nonlinear term 
G{u, dxu) and on the sign of the hrst eigenvalue of the linearized operator obtained from 
the linearization. Specifically, since we are taking into account also the nonlinear terms 
in V (arising from the linearization), we will show that both the form of G and the sign 
of Af influence the speed rate of convergence to zero of the perturbation. The estimates 
on V will be used to decouple the nonlinear system for the variables we end up 

with a one-dimensional equation of motion for the variable whose analysis is addressed 
in Proposition 3.11 In particular, the metastable behavior of the solution is described 


through the convergence of the interface location towards its equilibrium configuration. 
Again, the speed rate of this motion is influenced by the explicit form of G and by the 
sign of Af. 

Precisely, our results can be summarized in the following Theorem. 


Theorem 1.1. Let u{x,t) = U^{x;^{t)) + v{x,t) be the solution of the initial-houndary 
value problem (1.1) and let £ I such that U^{x‘,^*) is an exact steady state for (1.1). 
Then, for s sufficiently small, there exists a time , diverging to -|-oo for e —>■ 0, such 
that, for t < , the perturbation v is converging to zero with a speed rate depending on 
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£, and the interface location ^{t) satisfies the estimate 

|^(t) -ri < l^ole"^'* with /3= > 0 and /3" 


0 as e —>• 0. 


As a consequence, the interface location is converging towards its equilibrium configura¬ 
tion exponentially in time but, since is small in £, this convergence is extremely slow 
as e —7- 0. In particular, the solution u remains close to some non equilibrium configuration 
for a time that can be very long when e is small, before converging towards the steady 
state of the system, corresponding to 

This result characterizing the couple (^, v) gives a good qualitative explanation of the 
transition from the metastable state to the finale stable state. Also, since we are analyzing 
a complete system for the couple (^,u), the theory is more complete with respect to 
previous papers concerning metastability for parabolic problems (see, for example, |15L 
where only an approximation of the system is taken into account). 

Finally, in Section 4, we study, as an example, the case of a quasilinear viscous scalar 
conservation law: in this case we give an explicit expression for the approximated family 
{U^}, that can be used to provide an asymptotic expression for the speed of convergence of 
the interface, showing that it is proportional to Subsequently we analyze spectral 

properties of the linear operator arising from the linearization around the approximate 
steady state U^, proving that the first eigenvalue is negative and exponentially small in e 
(precisely of order while the rest of the spectrum is bounded away from zero. This 


analysis is needed to give evidence of the validity of the assumptions of Theorems 3.2 
3.41 and 13.71 at least in one concrete situation. 


3.3 


The main difference with respect to previous papers describing metastability for equa¬ 


tions of the form (1.1), and in particular with the work m, is that here we are considering 
a larger class of equations, where the form of the forcing term is not explicitly given and 
the diffusion is quasi-linear. The study of such class of equations could be a first step to 
address the problem of metastability for nonlinear-diffusion problems, such as the cases of 
the p-laplacian or the fractional laplacian diffusion operator. 

Moreover, as already stressed, in this paper we describe the behavior of the complete 
system for the couple (^,u), where also the nonlinear terms arising from the linearization 
are taken ion account. Since the forcing term G may even depend on the space derivative 
of the solution, we need an estimate also for the norm of the perturbation v. This 
gives a more clear overview with respect to m, since the complete system better suites 
the behavior of the solutions to (1.1). 


2. General framework and linearization 
Let us define the nonlinear differential operator 

V^[u] := £d^{a{x)d^u) - G{u,d^u), 

that depends singularly on the parameter e, meaning that V^[u] is of lower order. In 


particular, the evolutive equation (1.1) can be rewritten as 
(2.1) dtu = V^[u], u|^^Q = uo. 

Let us suppose that there exists at least one solution to V^^a] = 0, i.e. there exists a 
steady state for the problem (2.1), called here U^{x). Our primarily assumption is the 
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following: we suppose that there exists a one-parameter family of functions 0}?e/ 

such that 

\m-),v^[u%-,o])\ < E c{i), VC E I. 

being a family of smooth positive functions that converge to zero as e —)• 0, uniformly 
with respect to C- Moreover, we require that there exists a value C* E / such that U^{x; C*) 
is the exact steady state of the problem. 


The family can be seen as a family of approximate steady states 

for (2.1), in the sense that each element satisfies the stationary equation up to an error 
and that is measured by In particular, the parameter C describes 
corresponding to the location of the interface; if we suppose such 


that is small in e, 
the unique zero of 

parameter to depend on time, then the evolution of ^{t) describes the evolution of the 
solution to (2.1) towards its equilibrium configuration. 


Once the one-parameter family {U^{x',^)}^^i is chosen, we look for a solution to (2.1) 
in the form 


( 2 . 2 ) 


u{x, t) = U^{x-,^{t)) -I- v{x, t), 


where the perturbation v E {!)) is determined by the difference between the 

solution u and an element of the family of approximate steady states. By substituting 


(2.2) into (2.1), we obtain 


(2.3) 




dtv = Cy)V + r^[U%-,()] - — + Q^[u,C], 


where 

qv:=dr^[u^-,C)]v 

is the linearized operator arising from the linearization around , while Q^[u,C] collects 
the quadratic terms in v arising from the linearization and it is defined as 


Q^[u,C] :=V^[U^{-,0 + v]-'P^[U%-;0]-dVW{-,0]v- 


Example 2.1. Let us consider the case of a quasilinear scalar conservation law, i.e. prob¬ 
lem (1.1) with G{u,dxu) = dxf{u). We have 

C^'^v := edx {a{x)dxv) - dx := -\dx ■ 


On the contrary, when considering problem (1.1) with G{u,dxu) = g{u), since the forcing 
term G depends only on u, we obtain 

:= edx (a(x)dxv) - g'{U^)v, := -]^g''{U^)v^. 

In particular one has 


|Q^’^b,^]L < G\v\ 


Hi ■ 


< C\v\ 


l2‘ 


The form of the nonlinear terms in v will play a crucial role in the asymptotic behavior of 
the solution, as we will see in details later on in the is paper; in particular, it effects the 
speed rate of convergence of the solutions towards the asymptotic limit. 
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2.1. Spectral hypotheses and the projection method. We begin by analyzing the 
spectrum of the linearized operator we assume such spectrum to be composed of a 
decreasing sequence {^|(0}fceN of real eigenvalues such that 

• ^i(0 0 s-s e —0, uniformly with respect to 

• For all A: > 2, Z|(^) are negative and there holds 

Hence, we assume that there is a spectral gap between the first and the second eigenvalue 
and we ask for l\ to be small in e, uniformly with respect to 


Remark 2.2. We note that there are no requests on the sign of the first eigenvalue /f. 
Indeed, the metastable behavior is a consequence only of the smallness, with respect to e, 
of the absolute value of such first eigenvalue (see, for example, [M])- 

Denoting by ip\ = right eigenfunctions of and by ipf. = lAKsO the 

eigenfunctions of the corresponding adjoint operator we set 

Vk = Vk{^;t) := 

We now use an adapted version of the projection method in order to obtain an equation of 
motion for the parameter Since we have supposed the hrst eigenvalue of the linearized 
operator to be small in e, in order to remove the singular part of the operator in the 
limit e —7- 0, we impose ui = 0. Hence, the equation for the parameter ^(t) is chosen 
in such a way that the unique growing terms in the perturbation v are canceled out. In 
formulas 

and =0. 


|w(-;fW),A,t)) = o 


Using equation (2.3), we have 

+ S'lo.JI) + ■)§.»> = 0. 

Since, for small e, = Af(V'f,L’) = 0, we obtain a scalar nonlinear differential 

equation for the variable that is 

dt 


(2.4) 


We notice that if is the exact stationary solution, then 

vwi-, 0] = 0] - v^[u%--n] ^ qd^u%-,em - r )• 

Hence, at least for small £, the first eigenfunction ^|:f is not transversal to and we 

can take advantage of the renormalization 

(V>!(-;0,9gU^(-;0) = l, V^GI. 

Since we consider a small perturbation, in the regime u —)■ 0 we have 

i-Oq.o.v) = ‘ 

where the remainder R is of order o(|u|), and it is defined as 


R[v] := 


1 - 
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Inserting in (2.4), we end up with the following nonlinear equation for ^ 
(2.5) 
where 




p^[C,v] := ++ Q^[v,C])R[v]. 


Moreover, plugging (2.5) into ( |2.3[ ), we obtain a partial differential equation for the per¬ 
turbation V 

dtv = (£| + Ml)v+ TV^[v,^], 

where 


M\v 




3. The metastable dynamics 


The couple (u,^) solves the system 

r di 

(3.1) I dt 

{ dtv = + {Cl + Ml)v + n^[v, i], 

with initial conditions given by 


= 6*^(0 (1 + 


(^i(-;6),^^o - U{-;Co)) = 0, no = no - U{--,^o)- 


Our aim is to describe the behavior of the solution to (3.1) in the regime of small s. 


As stated in the introduction, the asymptotic behavior of the solution and, in particular, 
the speed rate of convergence of the shock layer towards the equilibrium configuration, is 
strictly related to the specific form of the nonlinear terms arising from the linearization 
of the original problem around the family U^. To be more precise, for a certain class of 
parabolic equations (as, for example, viscous conservation laws), these quadratic terms 
involve a dependence on the space derivative of the solution, so that an additional bound 
for the L^-norm of the space derivative of the solution is needed. On the contrary, when 
considering equations where the forcing term depends only on the solution itself (as, for 
instance, equation of reaction-diffusion type), we need to establish an upper bound only 
for the L^-norm of v. 

Furthermore, an important role is played by the first eigenvalue of the linearized op¬ 
erator; indeed, heuristically, the large time behavior of solutions is described by terms of 
order e^i*. In particular, the sign of Af characterizes the stability properties of the steady 
state around which we are linearizing. When such eigenvalue is negative, the steady state 
is stable and the solution is metastable in the sense that, starting from an initial config¬ 
uration located far from the equilibrium, the time-dependent solution starts to drifts in 
an exponentially long time towards the asymptotic limit. On the other side, when Af is 
positive, the solutions is said to be metastable because, starting from an initial datum 
located near the unstable steady state, the solution drifts apart towards one of the stable 
equilibrium configurations of the system, and this motion is extremely slow. 
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Hence, we need to distinguish different situations, depending on the type of equation 
we are dealing with; precisely, we will obtain different estimates for the perturbation v 
and for the speed rate of convergence of the shock layer (dictated by the behavior of 
solutions to (3.1), depending on the sign of Af and on the form of the nonlinear term Q^. 
Before state our results, let us recall the hypotheses we need. 


HI. The family {H''(-,^)} is such that there exist smooth functions such that 

(3.2) iv^l e cii), 

with converging to zero as e —)• 0, uniformly with respect to ^ G /. Moreover, we 
require that there exists a value G I such that the element U^{x;^*) corresponds to an 
exact steady state for the original equation. 


H2. The eigenvalues {A|(^)}j,gj, of the linearized operator are real and such that 
= .^i(0-4(6>ci and A|(0 <-^ for A; > 2. 

for some constants ci, C 2 > 0 independent on A: G N, e > 0 and ^ G /, and for some a > 0. 


H3. There exists a constant C > 0 such that 


|H^(e)l ^CIAKOI, vee/. 


3.1. The case Af < 0 and the quadratic term depending only on v. At first we 
consider the case of a nonlinearity that only depends on the perturbation u, and not 
on its space derivatives; we show that, if we consider a perturbation v such that u(0,x) 
is bounded, than we can perform an L°° estimate for the solution. In order to state an 
prove our result, we need an additional hypothesis. 


H4. Concerning the solution z to the linear problem dtz = C^z, we require that there 
exists > 0 such that, for all ^ G /, there exist a constant C such that 

\zit)\^,<C\zo\^,e-^^\ V^GI 

Remark 3.1. The constant C could depend on In this specific case, since ^ belongs to 
a bounded interval of the real line, if we suppose that ^ i—>• C'^(i) is a continuous function, 
then there exists a maximum of in I, namely C. 


Theorem 3.2. Let hypotheses Hl-2-3-4 he satisfied and let |uo|^oo < Too. Then, for 
e sufficiently small, there exists a time diverging to +oo as e —>■ 0, such that, for all 
t <T^ the solution v to 

\v\^oo{t) < C\Ll^\^,^t + 
for some positive constant C and 


(3.1) satisfies 


:= sup{Af(0} - 


> 0. 


The proof of Theorem 3.2 we present here is based on the theory of stable families 
of generators, firstly developed by Pazy in IIZI; it is a generalization of the theory of 
semigroups for evolution systems of the form dtu = Lu, when the linear operator L 
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depends on time. We refer to Appendix A for the definitions of the tools we shall use in 
the following. 


Proof of Theorem \3.^ First of all, let us notice that is a bounded operator that sat¬ 
isfies the estimate 




< ci|6»'^(0l < 


VC G I. 


(3.3) ||Af^ 

and is such that 

(3.4) 

for some positive constants ci and C 2 . Moreover, concerning the nonlinear terms and 
TZ^ and because of the specific form of Q^, there follows 

( 3 . 5 ) m,^<c\v\i^. 

Next, we want to show that is the inhnitesimal generator of a Cq semigroup 

7l{t, s). To this aim, concerning the eigenvalues of the linear operator we know that 
Af (C) is negative and goes to zero as e —)• 0, for all f, G I. Hence, defining Af := sup^ -^i(C)) 
we have A| < —|A|| < 0 for all A: > 1. Also, for t G [0,7"], is the infinitesimal 
generator of a Co semigroup 5^p)(s), s > 0 and, since H4 holds with the choice = |Af|, 
we get 

||%)(s)||<Ce-l^i^ 

and this estimate is independent on t. Thus, by using Definition 5.1 and the following 
remark, we can state that the family is stable with stability constants M = C 

and uj = — |Af|. Furthermore, since (3.3) holds. Theorem |5.2| states that the family 

T(i) 


In particular, by 


is stable with M = C and a; = — |Af | + 
choosing C = 1/C , oj is negative since H3 holds. 

Going further, in order to apply Theorem 5.3 we need to check that the domain of 
does not depend on time; this is true since -I- depends on time through 
the function C/^(x;C(A)), that does not appear in the higher order terms of the operator. 
More precisely, the principal part of the operator does not depend on 

Hence, we can define 7g(f, s) as the evolution system of dtv = (£| + AI|)u, so that 

(3.6) v{t) = 7^(f, s)uo+ f T^{t,r)H^{x] ^{r))dr+ f T^{t,r)Tl^[f^{r);v{r)]dr, 0 < s <t 

J S J S 

Moreover, there holds 

||rg(t,s)||<Ge-'^^(*-^), 


:= lA^I - C\Q^\^^ > 0. 


Finally, from the representation formula (3.6) with s = 0 and since (3.5) holds, there 
follows 


.(A) < e 


— I 


^o| 


, +sup \H^ 
C6L 


>(0 / e 

Jo 




dr+\v 


,{t) / dr, Vf>0, 

Jo 




t. 


so that, by using (3.4), we end up with 

{t) < e~^^^\vo\L^ + ci|D‘ 

Hence, setting N(t) := |r|^oo(A), we can rewrite the previous inequality as 

N{t) < AN^{t)+B, 
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and we can conclude N{t) < B providing AAB < 1. This condition is a condition on the 
final time t that reads 

(3.7) C,te-^^^^\vo\l^+C2mt^ <1. 

Precisely, the function g{t) := Cit behaves like |n^| for large t; 

since |n^| —>■ 0 as e —>■ 0, condition (3.7) is satisfied for all t < , where —)• oo as 

e —)• 0. 

Under this condition, the final estimate for v reads 

(3.8) |t'lioo(t) < +ci|U'^licx>t, for alH < 

and the proof is completed. □ 


Let us stress that, in the proof of Theorem 3.2, the negativity of Af is crucial in the 
construction of a stable family of generators. We also take advantage of the expression of 
Q^, where the first order space derivative of v does not appear; indeed, this allows us to 
estimate the nonlinear term TZ^ via the L°°- norm of v. 

If we start from an initial datum vq with a weaker regularity, precisely belonging to 
we can prove an estimate analogous to ( |3.8[ ) for the norm of the solution, but the 
following additional technical hypothesis is needed. 

H4.1 Given ^ G /, let and '!/’!(•;C) be a sequence of eigenfunction for the 

operators and respectively; we assume 

(3.9) 

j j 

for all k and for some constant C independent on the parameter 


Theorem 3.3. Let the couple {^,v) be the solution to initial-value problem (3.1). If the 
hypotheses Hl-2-3-4.1 are satisfied, then, for every s sufficiently small there exists a time 
such that, for every vq G L^{I) and for every t , there holds for the solution v 


V - z\ 2 {t) < C |U‘ 


+ exp 


A!(?(r))dr |uo| 


T2 


where the function z is defined as 

z{x,t) :=^ufc(0) exp(A|(^(r))(ir) ipl{x] f{t)). 


k>2 


Moreover, the time is of order \ sup Af(^)| ^, hence diverging to +oo as e 

i&i 


0. 


Proof of Theorem 3.3. Setting 


'^{x,t) = '^vfit) ip){x,i{t)), 


we obtain an infinite-dimensional differential system for the coefficients Vj 


(3.10) 
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where, omitting the dependencies for shortness, 

di- 


Fi ,= H^ + Y^ vj {^ ve) 


F2-.= 




(ScV’f, 


e „\2 




= -MF 


The coefficients aj, bjk are given by 

Oj := bji := (SjV'i, V^l) d^^j- 


Convergence of the series is guaranteed by assumption (3.9). Now let us set 


.E;fc(s,t) := exp / A|(^(T))dT 


Note that, for 0 < s < t, there hold 

Ek{s,t) = and 0 < Ek{s,t) < where A| := sup A|(0- 

Ek{0,s) ^6/ 

From equalities ( 3.10| ) and and since there holds vi = 0, there follows 

Vk{t) = Vk{^)Ek{^,t) 

+ [ {{'>PhH^) - + '^{il;l,bje)vejvj^Ek{s,t)ds 


+ 


[{'tPhO!') - {'^h-^^)]Ek{s,t)ds 

for k > 2. Now let us introduce the function 

z{x,t) := Ek{0,t) 


k>2 


which satisfies the estimate \z \^2 — I"^01 From the representation formulas for the 

coefficients v^, there holds 


v-z \ < 

' 'l2 


^^\^l^{m,F)\ + m,G)\)Ek{s,t)ds 


Moreover, since 


|0"(OI and m, H^)\ < C {^ + m,d^U^)\} 
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for some constant C > 0 depending on the L°°—norm of there holds 


[ \ ('^h^i)\Ek{s,t)ds 


k>2 j 

+ ^\{9ii^k,V’j)\\vj\ + ^\{'^k,d^^j)\\vj\^\{dii^l,V^e)\ \ve\^Ek{s,t)^ 

j j t 

<Cj]( rf^^(C)(l + |r;|2Ji?fe(s,t)ds)'. 

k>2 

On the other side, concerning the nonlinear terms, there holds 


Y^U' m,F2)\Ek{s,t)ds)'' < 

k>2 

k>2 ' 7 




1 - {d^i2l,v) 




Moreover, since < C\v\‘^^^, we have 


m,Qi\<c\v\ 


iw.e‘>i , ,on-.|2 


|X 17/1 J- 

\e^^^rl,vf{rk,^^u^)\<cn%o\vt 


so that we end up with 


Y{fMk,F2)\E,,{s,t)dsy<CY{fK,i'^ + ^"(0) Ek{s, t) ds 

k>2 do 7.>o -^0 


Since y/a + b < y/a + Vb, we infer 


-IZ / + + 1^1%} Fk{s,t)ds 

k>2 

<C f |o"(0(l + |p|JJ + }E Ek{s,t) ds. 
d^ k>2. 
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The assumption on the asymptotic behavior of the eigenvalues A| can now be used to 
bound the series. Indeed, there holds 

Ekis, t) < E2 {s, t) Y ^2(s, t). 


k >2 


k >2 


As a consequence, we infer 


u — z| 2 < C* / VE{^){t — s) E2{s,t) ds 


+ C y ||u - z|^2 + I^li2} E2{s,t) ds. 


Now, setting N{t) := sup iiii(s,0) |(u — ^)(s)| 2 ) we obtain 
se[o,t] ^ 

Ei{t, 0 )\v-z\. <C [ N^{s){t-E 2 {s,t)Ei{ 0 ,s)ds 


+ C E2{s,t) Ei{s,Q) ds 

Jo 

+ C j {t-E 2 {s,t)Ei{s,Q)ds 

< CiN\t)E^{Q,t) + + |A||-3/>o| 


2 

l2 ; ’ 


where we used 


(3.11) 




{t-s) E2{s,t)ds< [ {t-s) V2gA|(t 


JO JO 

and Cl and C 2 depend on A^. Hence, as soon as 

(3.12) 4Ci (|A||-V2|q^|^^ +s,(0,t)|A||-3/>oi;2) < 1 
we obtain the following L2—estimate for the difference v — z 

(3.13) \v - 2|^2 < |uo|^2 '^i(OT)) , 

where Ei{0,t) behaves like since Af < 0, and [All ~ e““ for some a > 0. Condition 
(3.12) is a condition on the final time and it can be rewritten as 


< C 


1 - 


|A||-3/4|uo|^2 ' 


Hence, can be chosen of order (inlAlP^^) |A^| which is diverging to +00 as |Af| ^ 
for e —)■ 0. 

□ 
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3.2. The case Af > 0 and the quadratic term depending only on v. Since Af 
is positive, we can no longer use the theory of m-, indeed, we cannot state anymore that 
< — |A^| < 0 for all fe > 1, so that we cannot construct a stable family of generators 
for as in the proof of Theorem 3.2 In this case we can prove an estimate analogous to 


the one proved in Theorem 3.3 


Theorem 3.4. Let the hypotheses of Theorem 3.3 hold, 
holds for the solution v 


Then, for every t < T^, there 


(3.14) 


\v — z 


l2 




where the funetion z and the time T^ are defined as in Theorem 


3.3 


and 1^2 :=supA 2 ( 0 - 
i&i 


Remark 3.5. From hypothesis H2, we know that Al < for some a > 0. Hence, 

if OL is strictly positive, (3.14) assures the convergence to zero of the perturbation v as 
£ —0; on the contrary, if a = 0, we need to restrict our analysis to the case of small (with 
respect to e) initial data, i.e. we need to require vq G L‘^{I) such that |uo |^2 ^ ce. This 
is a small deterioration of the estimate, consequence of the instability of the steady state, 
but it is however consistent with the case considered in |6]. 


Proof. The proof is exactly the same as in Theorem 3.3; recalling (3.11), we end up with 
the following estimate 


\v — z\ o < 

I Ii2 — 


(|n'|,„+|A|r=>/-|t,„|2^Ei(0,*)). 


Since A^ is positive, we can no longer assure the convergence to zero of the term |?io |^2 

0, Ei(0,t) is converging to a constant for small 


0 as e 


in (3.13); indeed, since Af 
s and large t. We here use 
recall that, if a = 0, we have the assumption |'yo |^2 < ce. 


e and large t. We here use hypothesis H2 concerning the behavior of A|, /c > 2 and we 


□ 


Remark 3.6. Let un consider the Allen-Cahn equation, i.e. problem (1.1) with a{x) = 1 
and G{u,dxu) = g{u), for some g satisfying the following assumptions; there exists a 
potential function VF : M —>■ M such that g{u) = W'{u), and we assume W to have two 
distinct global minima Pu* such that W{P.u*) = 0; in this case, is its well known (see, 
among others, mm) that the only possible stable equilibrium solutions are constant 
in space and are given exactly by ±u*, while all the space dependent steady states that 
present patterns of internal transition layers are unstable. Additionally, in [6], it is proven 
that the first eigenvalue of the linearized operator around an interface stationary solutions 
(which is unstable) is positive, but small in e. Also, the nonlinear terms depend only 
on the perturbation v, and can be estimated via the L^- norm of v. This is an explicit 
and well known example where the equation exhibit a metastable behavior and Theorem 


3.4 can be applied (see, for instance, 

Our guess is that, also in the case of a quasilinear second order term with a{x) satisfying 


(1.2), a spectral analysis of the linearized operator can be performed in order to show that 


Af is positive and small in e. 
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3.3. The case depending on v and its space derivatives. 

In order to prove an estimate for the perturbation v, we need an additional upper bound 
for the L^-norm of dxV and we have to consider initial data vq. 


Theorem 3.7. Let hypotheses H1-2-H4.1 he satisfied and let us denote by {^,v) the 
solution to the initial-value problem (3.1), with 

C( 0 ) = Co & I and 


u(x,0) = uo(x) G 

Then, for e sufficiently small, there exists a time > 0, such that, for any t < , the 

solution V can be represented as 

V = z + R, 

where z is defined by 


z{x,t) := y^Vk{0) exp 


k>2 


and the remainder R satisfies the estimate 


(3.15) 


I^Li <c 


exp 


Ai(C(F))dr |uoL. +e 


1-5 




for some constant C > 0 and for some 6 G (0, a), a > 0. Furthermore, the final time 
can be chosen of order 1/e^, for some 7 > 0 . 


Remark 3.8. As will be clarified later, the constant a in (3.15) is exactly the one defined in 
hypothesis H2. We do not consider the case a = 0 since it appears only when considering 
reaction diffusion systems, where the nonlinear term depends only on the function v 
(see, for example, [23]). 


Remark 3.9. Since a > 0, the final estimate for v in Theorem |3.7| does not depend on 


the sign of the first eigenvalue Af; indeed, the term e^£'i(0, t)|uo|^^ goes to zero as e —)■ 0 
whatever the sign of is. Hence, with respect to Theorem |3.4[ Theorem 3.7 holds for a 
general class of initial data vq G H^, not only the ones with norm small in e. On the 
contrary, we also underline that the estimate (3.15) for v is weaker that the corresponding 
one obtained in (3.13) in Theorem 3.3 indeed, it states that the remainder R tends to 
0 as instead of | | , and we shall see in the following that this term behaves like 

Such deterioration is a consequence of the necessity of estimating also the first 
order derivative. 


Proof of Theorem 3.1. Since the plan of the proof closely resemble the one used for proving 
Theorem |3.3[ we propose here only the major modifications of the argument. In particular, 
the key point is how to handle the nonlinear terms. 

Setting as usual 

= Y^Vj{t)(p‘jix,C{t)), 

j 

we obtain an infinite-dimensional differential system for the coefficients Vj 


dvk 

dt 


^kiO '^k + (V’fc) -^1) + (V’l) -^2) 
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where Fi is defined as before as 

Fi:=H^ (“i + Yj '^e) 

j e 

with 

aj := , (f^.) bje := p>\) 

The term F 2 comes out from the higher order terms and TiF and has the following 
expression 






,s ^\2 




F2 := Q" - I - 

Moreover, we have 

l(i/'|, 0 )| < (1 + 

so that, setting as before 

Ek{s,t) :=exp^^ A|(^(r))dr^ 

and since vi = 0, we have the following expression for the coefficients Vk, k > 2 

Vk{t) = Vk{0)Ek{0,t) + f {{i;l,Fi) + {'il;l,F2)}Ek{s,t)ds. 

Jo 

By introducing the function 

z{x,t) := '^VkiO) Ek{0,t) pl{x;C{t)), 

k>2 

we end up with the following estimate for the L^-norm of the difference v — z 

rt 


V - z 


-Y [ + + Ek{s,t)ds, 

k>2 


that is 

\v — z\ „ FC f — s)~^^‘^E2{s,t)ds 

Jo 

+ Jo - ^1^1 + l^l^i) - s)-^^^E2{s,t)ds, 

where we used 

Y,Ek{s,t)<C{t-s)-^/^E2{s,t). 

k>2 

Now we need to differentiate with respect to x the equation for v in order to obtain an 
estimate for \dx{v — z)\^ 2 - By setting y = dxV^ we obtain 

dty = C^y + + dx {dV^U^]) v + dxH^x,^) + dxF^[v,^], 

where 


:= —dx^U^{-;^) 6*^(0 dx {dV^[U^]) v := ea'{x) dxV — dG[U^]v. 
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Hence, by setting as usual 




we have 


where 


^ = m) vk + m. F*) + d, {dr^[u^]v)) + m, d.n^), 




1 + ^ viid^tpl, Pi) 




Moreover, by integrating by parts, for some m > 0 there hold 


Wk,d.n^)\<C\v\l^, (V-|,5.(dP^[C/>)) +e|a'|2^) + — 

Because of the assumptions on a and G, there holds + |a^|^oo < c for some positive 

constant, so that, by integrating in time and by summing on k, we end up with 

\y-d^z\^^ <C |f7^(0(l + |w|j2)++ +e™| Ei(s,t)ds 

+ F j l^l^i + X/ 


Now, given n > 0, let us set 


1 


F^{t) ■=-p sup \v - Ei{s,0), 

^ se[o,t] 


so that we have 


(3.16) 


and 

(3.17) 


^Ei{t,0)\v - z\^ < C f ^^-^(t-s) ^/^E2{s,t)Ei{s,0)ds 
£ Jo ^ 


+ ^ +l^l^i) (i-s) ^^^E2{s,t)Ei{s,0)ds. 


^Ei{t, 0 )\y-d:^z \^2 <C j |£;i(s,0) + (t - s) ^^^Esis,t) Ei{s,0)'^ ds 

11 

Ei{s, 0)ds 


+ C 


+ c 


1 1 

— + 


u - zP + IzP I + 


1 1 

— + 


^n+m 




J/I ' 'j/i/ 


^|(t-s) ^^‘^Es{s,t)Ei{s,Q)ds. 
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By summing (3.16) and (3.17) and since there hold 

1 


rt 


rt 


JO Jo 

we end up with the estimate N{t) < AN‘^{t) + B, with 


|A||i/2^ 


+ e-^-"^|A||->o|^i + £^-^Ei{t, 0) {t + |A||-V2). 

We now use the assumptions on the behavior of A|, /c > 2; since there holds jAH ~ e 
for some a > 0, if we require m < a, for all n > 0 there holds N{t) < B that is 


(3.18) \v-z\^,<Cm^^ + [e'--^\vo\l^E,{0,t)+e^). 

Precisely, we can choose m = a — 5, for some S G (0,a). Finally, providing m > n, we 
can choose the final time of order 0{e~'^) for some 0 < 7 < 1. Now the proof is 

completed. □ 


Remark 3.10. The proofs of Theorems 3.3 and |3.7| can be easily extend to the case 
V G [77^(1)]"' with only minor changes; this is meaningful in light of a possible application 
of these results in the case of a system of equations of the form (1.1). We point out that, in 
the case u G [77^ (/)]"■, in hypothesis H2 one should consider the chance of having complex 
eigenvalues; however, the first eigenvalue has to be real in order to observe a metastable 
behavior. 


3.4. The slow motion of the shock layer. The estimates for the perturbation v ob¬ 
tained in the previous section can be used to decouple the system (3.1) in order to obtain 
an equation of motion for the parameter ^(t). Indeed, both (3.8), (3.13) and (3.18) im¬ 
ply that, for small e, the perturbation v converges to zero as 7 —?• 00 . This preludes the 
following result. 


Proposition 3.11. Let either the hypothesis of Theorem 3.2, 3-4 
us also assume that 


or 


3.1 be satisfied. Let 


(3.19) (^-r)^^(0<0 foranyiGlA^fi* and r'(r) < 0. 

Then, for e sufficiently small, the solution ^(7) converges exponentially to as t ^ -|-oo. 

Proof. Again, we need to divide the proof in two main parts. At first, we consider the 
case where the nonlinear higher order terms in the equation for the perturbation v can be 
estimated via the L^ norm of v itself. Since either ( |3.13 ) or ( |3.14 ) hold, we get 

§ = 9'(0(l+r)+p' 


with |r| < C -b |A|| and |/| < CIAH 
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By a method of separation of variables, and since 0^{^) r-^ — we end up with 

m ~ r + (Co - + c\m-^/^\vo\% (i - 


where /3^ := —+ |r|) ~ for small e, and corresponds to the asymptotic 

location for the parameter ^ (we recall that U^{x,^*) is an exact steady state for the 
system). Hence, for small e and large t, ^ is converging to with exponential rate. 


When the nonlinear term also depends on the space derivative of v, we use (3.18) and 




we get ^ = 0^(^)(l+r)+/9'^, where both r and can be estimated via (e”^+e 
Again, since both r and go to zero as e —)• 0, in the vanishing viscosity limit there holds 

showing the exponentially (slow) motion of the position of the interface towards its equi¬ 
librium location . □ 


4. Application to quasilinear viscous scalar conservation laws 

In this Section we mean to apply the general theory previously developed to the specific 
example of a quasilinear conservation law, i.e. we consider the following initial-boundary- 
value problem 


(4.1) 


X G /, t > 0, 


t > 0, 


dtu = edx {a{x)dxu) - dxf{u), 
u{±£, t) = u±, 

u(x, 0) = uo(x), X G/, 

where the flux function / G C^(I) satishes the standard hypotheses 
(4.2) f"(u)>co>0, f'(u+) < 0 < f'(u-), f(u+) = f(u-). 


This is a well know and simplified prototype of problem (1.1) in the case of a forcing term 


G(u, dxu) depending both on u and on its space derivative. 

Formally, in the vanishing viscosity limit e —>■ O"*", equation (4.1) reduces to a first-order 
quasi-linear hyperbolic problem on the form 

(4.3) dtu +dxf{u) = {), ri(x,0) = tto(x) 

complemented with boundary conditions 

(4.4) u{—i,t) = u- and u{i,t)=u-^-. 


The standard setting of solutions to (4.3) is well known, and it is the one given by the 


entropy formulation. Hence, we may have solutions with discontinuities, which propagate 
with a speed s dictated by the well known Rankine-Hugoniot relation 

sM = I/(^^)l, 


and that satisfy appropriate entropy conditions. Assumptions (4.2) on the flux function / 


guarantee that the jump from the value u- to the value u+ is admissible if and only if 
u- > u+, and that its speed of propagation s is equal to zero. 
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U- 


/ 

entropic jump 


\ 


+ 1 


non entiopidJump 


'u+ 


Figure 1. Picture of the steady state to (4.3); because of the entropy condi¬ 
tions, the only jumps admitted are the ones from a value U- > 


In this case, equation (4.3) admits a large class of stationary solutions satisfying the 
boundary conditions, given by all that piecewise constant functions in the form 


u{x) 


u- X e {—I, xo), 
u+ x£(xoJ), 


where xq G / is a certain point in the interval (see Figure 1). Hence, given ^ G we 

can construct a one-parameter family {b^i,yp('iO} of steady states, parametrized by ^ that 
represents the location of the jump, and given by 


C^hyp(^;0 = u-X(-e,o(^) + u+X{^,£){x), 


where xi denotes the characteristic function of the interval I. 

For the initial-boundary value problem (|4.3|)-(4.4), it is possible to prove that, starting 


from an initial datum uq with bounded variation, every entropy solution converges in finite 
time to an element of the family {t4yp(s0}- 


For e > 0, the situation is very different. In this case, because of the diffusive term, no 
discontinuities are admitted, and there is a drastic reduction of the number of stationary 
solutions; indeed, it is possible to prove (see Section 4.1) that there exists only one single 
steady state satisfying 

I edx {a{x)dxu) = dxfiu), 

I u{±i) = u±. 

Such solution, denoted here by U^{x), converges pointwise in the limit e —0+ to a specific 
element t/i,yp(-;C) of the family {t4yp(s0}) for some G I. 

Finally, the single steady state is asymptotically stable (for more details see the 
spectral analysis performed in Section 4.2), i.e. starting from an initial datum close to 
the equilibrium configuration, the time dependent solution approaches the steady state 
for t —)• -boo. 
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4.1. The stationary problem. In order to apply to problem (4.1) the general theory 
developed in Section 3, the first step is the construction of the family of approximate 
steady states defined in hypothesis HI. 

Let us then consider the stationary problem for ( |4.1| ), that is 


(4.5) 


e dx {a{x)d^u) = dxf{u) 
n(—£) = u_, u{i) = 


for X G {—£,£). This problem has been extensively studied in the case a(x) = 1 (see, for 
instance [IQ]). Let us begin our study with some explicit examples. 

Example 4.1. The Burgers equation. In the case of Burgers equation, i.e. f{u) = tt^/2 
and a{x) = 1, the value coincides with —tt_, and the stationary problem reads 

edlu = ]^dx{u^), u{±t) = =Fu*. 

By integrating, we obtain an explicit expression for the unique steady state, that is 

L’^(x) = —K tanh > 

where k = K.{e,i,u*) is univocally determined once the boundary conditions are imposed. 


Following the general approach introduced in the previous sections, we want to construct 
the family of approximate steady states {U‘^{x]^)}. There are several choices to built-up 
such a family (for instance, the approach of [7j based on the existence of traveling waves on 
the whole line). We here recall the approach of |15) . and we consider a function obtained 
by matching two different steady states satisfying, respectively, the left and the right 
boundary condition together with the request U^{^) = 0; in formulas. 


(4.6) 


tanh (k_(^ — x)/2e) in (—^, 0 

K+tanh (k+(^ — x)/2e) in 


where k± are chosen so that the boundary conditions are satisfied. By direct substitution 
we obtain the identity 

in the sense of distributions, where 6x=^ the usual Dirac’s delta distribution centered in 
X = ^. Going further, we have 

= ^(k- - k+){k- + K+). 

In order to determine the behavior of for small e, we need an asymptotic 

description of the values k±. Let us set k± := =F^±(1 + h±) and A± := ^ =F Then it 
results 

ta„h(T^(l + /i±))=^. 

Therefore, the values h± can be chosen both positive and then 


0 < tanh 



u±^± '\ 

2e ) 


1 

^l + h±' 


2e 
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that gives the asymptotic representation 

1 


h±<- 


- 1 = 


= 2e^“±^±/^ + l.o.t 


tanh (=Fri±A±/2e) g=F«±A±/£ _ ^ 

where l.o.t. denotes lower order terms. In particular, since u± = ^u* for some u* > 0, 
there holds 


h± 


^-u*A±Ie 


Finally, since 
we end up with 


1 

ldxU%=^ = —(fe- - k+){k- + A:+)~y (/i_ - h+), 




showing that this term is exponentially small for e —)• 0 and it is null when = 0, that 
corresponds to the equilibrium location of the shock when f{u) = v?I2. Hence we have 
the following asymptotic expression for the term defined in 






that shows that hypothesis HI is satisfied in this specific case. 

Example 4.2. The quasilinear viscous Burgers equation. In the case of a generic function 
a{x) satisfying (1.2), the expression for the unique steady state is given by 

(4.7) 

where 


U^{x) = —K tanh 


K [b{x) - 6(0/2] \ 
2e 


h{x) := J 


1 


-e a{t) 


dt, 


is known to exist because of the assumption Let us observe that, if we call N the 

number of zeroes of b{x) — b{l)/2, the number of the layers of the steady state U^{x) 
defined as in (4.7) is exactly N; hence, since the only jumps admitted are the ones from a 
a value ui to a value U 2 < ui, and since the boundary conditions are such that U- > tt+, 
there exists a unique solution to (4.5) if and only if the function b{x)—b{l)/2 has a single 
zero located at some point x* G I. In this case, the steady state will be approximately 
given by 

The family of approximate steady is constructed as in Example |4.l| recalling that now 
is null when x = x*. Hence 


U%x;0 = 

In particular, there holds 


tanh (k_ [(^ — X*) — x] /2e) 
K+ tanh (k+ [(^ — X*) — x] / 2e) 


in (-£,^-x*) 
in (^-x*,£), 


so that 
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which is null when ^ = x*, corresponding to the location of the interface of the exact 
steady state for the problem. 


In the general case, if the flux function f{x) satisfies hypotheses (4.2), solutions to (4.5) 
can be found implicitly via de formula 

r- ds 


(4.8) 


'i{x) « - f{s) £ J- 
where k G {f{u±), +oo) is such that 

ds 1 


1 j-X 

- / a-\x) 
£ J-£ 


dx 


4>(k) := 


iu+ 


- f{s) £ 


a ^(x)dx 


Assumptions (4.2) on the flux / imply that ‘h is strictly decreasing and such that 

lim 4 >(k) = + 00, lim 4 >(k) = 0 . 

K—K—>- + oo 


Therefore, for any > 0, there exists a unique solution to (4.8) satisfying the boundary 
conditions. 

As in the previous example, the family is build up by matching at ^—x* G {—i, 1) 
two different steady stated Ut and I/^, solutions to the stationary equations in (—^, £,—x*) 
and {^—x*,i) respectively; hence we have 


(4.9) 


and 


U%x;0 = 


Ul{x]^) —i<x<^ —x*<i 

U^{x; 0 —£<^ — x*<x<£, 




and - x*-,C) = u*, Ul{£;C)=u+, 


where u* is such that f'{u*) = 0 and x* is the unique zero of the exact steady state C7^(x). 
Moreover, without loss of generality, we assume that f{u*) = 0. Note that rt* = 0 in the 
case of a Burgers flux f{u) = v?/2; this is consistent with the interpretation of ^ as the 
location of the interface of u. Once again the error V^\U^] is given by 

Now Ul{x]^) and U^{x;^) are implicitly given by 

ds r _i 


Ju‘ 

f 


h rr ^ xG{-l,i-x*) 


ds 


a '^{s)ds xG{x — x*,1) 


- f{s) J^-x* 

where k± are chosen such that the boundary conditions are satisfied. Then 

^a-^{x)if{Ui{x-0)-k.) 


dxU%x]() = < 


a-\x){f{Ul{x-,0)-k, 


X G {—I, ^ — X*) 
X G(^- x*,l), 
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and we immediately obtain 




1 


ea(^ — X*) 


ifmiC - x*;0) -k+- - x*-i)) + k.) = 


= —T^ -^ (fc- - k+) < —\k- - k+\. 

£a[i-x*)^ - ae' 

Using the following bounds on / 

f{u±) + f{u+){u - n+) < f{u) < {u^ -u) ue [u+,u*], 

u* — M+ 

f{u±) - f'{u-){u- -u) < f{u) < - u*) uG [u^,u-] 

U- — u* 

we can proceed as in [15] to estimate |/c_ — k+\, proving that, for any 6 G (0,£), there 
exists C > 0, independent on s, such that 

£\ldxU^^ 

showing that hypothesis HI is satisfied. 

4.2. Spectral analysis. We mean to analyze the spectrum of the operator 

Cp := edx {a{x)dxv) - {f'{U^) v) , 


= ^ —a:* I — 


G {-I+ 5,1-6), 


obtained from the linearization of (1.1) around an element of the family (4.9); in particular, 
we obtain a precise distributions of the eigenvalues A|(^). 

This analysis is needed in order to show that the general theory previously developed 
is applicable in the specific case of quasilinear viscous scalar conservation laws; more 
precisely, we show that the hypothesis H2 concerning the distribution of the eigenvalues 
of the linearized operator is satisfied in a concrete situation. 

The eigenvalue problem reads 

£dx{a{x)dx+) - dx{f'{U’^)+) = X^(p, ¥^(0) = ip{£) = 0. 

Firstly, we show that the eigenvalues of are real. To this aim, let us introduce the 
self-adjoint operator 


where 

(4.10) 


:= £dx{a{x)dx'ip) - W^x] 


w%x-,m) ■■= 


1 

a{x) 


+ ^edxfm. 


A straightforward computation shows that (p^ is an eigenfunction for £| relative to the 
eigenvalue A*" if and only if 

= exp (-1 £ 

is an eigenfunction for the operator J\f^ relative to the eigenvalue = eX^. Hence 
(4.11) ea(£|) = a(AA|), 

so that, since is self-adjoint, we can state the the spectrum of is composed by real 
eigenvalues. 
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Going further, if u is an eigenfunction of relative to the first eigenvalue Af, integrating 
in the relation = Af u, we deduce the identity 

ri pi 

0 = / — X\)udx = e {a{l)u'(1) — a{—l)u'(—1)) — Xl / u{x)dx. 

J-e J-e 

Assuming, without loss of generality, u to be strictly positive in {—£,£) and since a{x) > 0 
by assumption, we get Af < 0. Hence, there holds 

cr(£|) C (-00,0). 

Remark 4.3. With analogous computations, it is possible to prove that the eigenvalues 
of the linearized operator obtained from a linearization around the exact steady state 
U^{x) are all negative; this shows the asymptotic stability of U^{x). 


4.2.1. Estimates for the first eigenvalue. We mean to control from below Af. To this aim, 
we estimate the first eigenvalue ^f of the operator and we use the relation (4.11); by 
means of the inequality 


\hl\< 


IV'I 


l2 


that holds for smooth test function such that fiiH) = 0, we look for a test function 
such that fi{x) := iPq{x) — K^{x), where 

^SW:=exp(l^ 

and such that there holds 

:= 

A direct computation, show that has to solve 

r dx {a{x)dxK^) = 0, 

1 K^±e) = roi±^). 


Hence, by integrating, we get 


K%x) := 


a{x) 


-1 


r dy 
Li a{y) 




Going further, there holds 

w< 


<c 


\K^\ 


C 


as soon as |'!/’ol^2 > l-^^l^2- We assume the opposite case being similar; 

from the definition of and from the properties of a(x), it follows 

<2 £i/; 2(-£). 


\K^\\ 

' 'l 2 


so that 


\K‘ 




e\-2\jR\2 
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We define 


r := 




Since converges to the step function := u- X{^,£) as e —)■ O"*", we get 


P = 


exp 


i-i 




e J-e a{y) 


L 


e J-i 


nu\i\y)) 

a{y) 


dy ) dx 




> e 


where we used a{x) < /3. Moreover there holds 


/0 = 


exp 


fi\ fin-) ^ r !M 1 A\ 

Ve yj-e a{y) ^ a{y) ^\) 


dx 


> 


ft 

/ e^' 

J-e 


(u ){i+t)/p£ +f'{u+){x-^)//Se 


J-t 


s. 


In particular, if we suppose \f{U^) — f'{U^)\ i < cqs for some cq > 0, we end up with 


\K' 


e| —1 l„l,£ 




1^2 IY0li2 


Thus, for the first eigenvalue y\ of the self-adjoint operator there holds the estimate 


> ( 
idj' 

\yW < [Ci y/e for some positive constant Ci,C 2 - As a consequence, since the 

spectrum <t(£|) coincides with e“^cT(A/’|), there holds 

(4.12) 

Remark 4.4. Since /' is a continuous function, the request |/'([/^) — < cqE is 

satisfied if we require \U^ — < cqE, which is consistent with the convergence of 

to as e ^ 0. 

4.2.2. Estimate from above for the second eigenvalue. We mean to give an estimate on the 
behavior of the second and subsequent eigenvalues of the operator To this aim, we need 
some additional assumptions on the limiting behavior of the functions a^{x) =f'{U^{x)) as 
e —>■ O'*'. Precisely, inspired by [15], we suppose that G ([—£,£]) and a G 
satisfy the following hypotheses: 


i. g CH[-£,£] \ {e}), a G C\[-£,i] \ {C}) and 

da^ dfa^ ^ da , „ 

— ,^<0<a",— m 

dx dx^ dx 


da^ da 


2„£ 




in (^,£), 


' dx' dx dx‘^ 

ii. there exist the left/right first order derivatives of a^ at ^ and 
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da^ 


dx 


< C 


and 


liminf e 

£— 




> 0 . 


iii. for any C > 0 there exists Cq > 0 such that, if |x — ^| > c^e, then 

|/'([/^)-/'([/°)| <Ce, 

where we recall f{U^){x) := + f{u+)x{^,e)- 

Under these assumptions, it is possible to prove the following Lemma describing the 
function + eX^, with given in (4.10). 


Lemma 4.5. Let the family and the function a he such that assumptions (1.2), Al-2-3 
are satisfied, and let < 0 6e such that 

inf eA^ > where := min{|/'(u“)|, |/'(u+)|}. 

Then there exist > 0 such that, for e < e^, the function + eX^ enjoys the following 
properties: 

i. + eA^ is decreasing in and increasing in 

ii. there exist C, c > 0 such that, for any x with |a: —> ce there holds +£X^ > C > 0; 

iii. there exist the left/right limits ofW^ + eX^ at ^ and ft := limsup(lU^(^±) +eA^) < 0. 

£—>0 + 

Proof. The proof lies on a straightforward application of the properties of the functions 
f'{U^) and o(a:), and closely resemble the one of [l5l Lemma 4.3]. □ 

Remark 4.6. In the easiest case of the Burgers equation, i.e. a(x) = 1 and f{u) = u‘^12, 
from the explicit expression of given in (4.6) and since f'{u) = u, it is easy to check 
that the assumption we made on f'{U^) are satisfied (see also Figure 2). 


From Lemma 4.5 we can infer that, in the regime e —)• O'*', + eA*" has two zeros in 

[—f,£], denoted here by moreover 

-(■ < yt < i < y% < I, and |y| - ^| < e. 

Let Al and pL 2 = sX^ be the second eigenvalues of the operators and A4| respectively, 
with corresponding eigenfunctions y?! and 'i/’f such that 

(4.13) =exp(-l J). 

and let assume that Lemma 4.5 holds with A^ = A^- We stress that we can choose A^ = A^ 
without loss of generality since, if not possible, then it would follow A^ ~ —l/e“ with 
a > 1, which implies that H2 is trivially satisfied. 

Since Al is the second eigenvalue, applying the Sturm-Liouville theory to the operator 
A/"!, we deduce that the functions (/?! and '01 possess a single root located at some point 
Xq G (—£,£). The sign properties of + //I described in Lemma 4.5 imply that Xq G 
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Figure 2. The approximate steady state for the Burgers equation. C/® is ob¬ 
tained by matching two exact steady states in the intervals and (^,^); as a 

consequence, t/® is a function but its first order derivative has a jump located 
in a; = ^. 

(yi,y+). Indeed, by contradiction, it is easy to verify that if Xq G yi) U then 

it would exists at least one point x G yl) U {y%,£) such that 

dxil’2ix) = 0, / 0, V'K®) 7^ 0, 

and, using the equation solved by this would imply that d'^'ijj^ix) and V’ 2 (^) have the 
same sign, which is not possible. Now (^1 V’l restricted to the intervals {—i, Xq) and 
{xQ,i) are eigenfunctions relative to the first eigenvalue of the same operator considered 
in the corresponding intervals and with Dirichlet boundary conditions. 

Without loss of generality, we can assume Xg > ^ and y?! > 0 in (xq,£) and we can 
restrict our attention to the interval J = (xg, t). By proceeding as in [15] , integrating on 
J, we get 


fi 

A| / (t)ldx = e {a{l)d^(j)l{t) - < -ea{xl)d^(l)\{xl). 

Jx% 


If we now assume '01 lo be as in (4.13) with xg = Xg and renormalized so that max 01 = 1) 
from the previous equality we deduce 


(4.14) 

where 


Is ■ = 


IA 2 I > ea(x^) 5a;0|(x^) 4 0 

nu^iy)) 


exp 


0 


a(.y) 


dy dx. 
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In order to get an estimate from below on | |, we give an estimate from above on and 

an estimate from below on 


le < e 


< 


|a"-a°|^l/2ea {u-){x-xl)/2 sol ^ l^i/2£a 

Jxf, l/'(^“)l 


2ea Ja‘"-a°|^^/2£a ^ ^ ^ 


\f{u )| 

where the last inequality holds since la*" — a*^| 


<Ce. Hence, (4.14) becomes 


(4.15) 


|A||>cf(4) 


for some C > 0 independent on e. Now, let xm £ be such that 'ijj^ixM) = 1; from 

the properties of + e X stated in Lemma 4.5, it follows that xm £ ixo,y+)- Then, by 
Lagrange Theorem, there exists xl G {xq,xm) such that 


dV’2 




1 

-> 

XM - xl 


1 


> 


1 


dx XM - xl y+-C 

Since iplixl) = '02(^) = 0 > 0 ™ (^ 0)^)5 infer that the function tpl is 

concave in the interval (xQ,y+), deducing that 


f ("»> 2 


> - 
dx 


CnE 


In conclusion, from (4.15), we deduce 
(4.16) |A|| > 


C 


A|<- 


C 


for some C independent on e. 

Estimates (4.12) and (4.16) show that hypotheses H2-H3-H4 are satisfied in the case 
of a quasilinear viscous conservation law. 

4.3. The speed rate of convergence of the shock layer. We here mean to obtain an 
asymptotic expression for the term 0^(^); indeed, recalling the equation for ^ in (3.1), and 
supposing the perturbation v to be small, we have 

f 

Hence, the function 6^ gives a good approximation of the speed rate of convergence of the 
solution towards its asymptotic steady state. 

Without loss of generality, we may assume x*, the unique zero of the exact steady state, 
to be equal to zero. 

Since 6^{^) := (V'f (0) ^"^[17^]), we need an expression for i/;f, the first eigenfunction of 
the adjoint linearized operator defined as 

:= edx {a{x)dxv) + f'{U^)dxV. 

Following the idea of m. for e ~ 0, the eigenfunction 'ijjf is close to the eigenfunction of 
C^’* relative to the eigenvalue A = 0, where 

/'([/°)(x; 0 := (a:) + /'(m+)x., ,, (x) 
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Hence, V’i(®) solves 

edx{a{x)d^'il;i) + = 0, for x G 

edx{a{x)d^'il;°) + f{u+)d^'il;° = 0, for x G 

.V'?(-^) = 0, '0?(^) = O, = 0. 

By integrating in and respectively, and by imposing the boundary conditions 

and the condition on the jump, we obtain the following expression for 


^\{x) := 


(1 

(1 


_ p-f'{’»■-)W)+Ki))h^{i _ pf'{'^+)W)-K^))h\ 


x<i, 

X> 


being b{x) ■= f a ^{x) dx. In particular, in the limit e —)• 0, we obtain ss 1 so that 

This estimate show that the speed of the interface is exponentially small when e is small; 
for example, in the special case of the Burgers equation there holds 




^ g-«*(£+0/£ _ g-“.(^-?)/£ 


showing that the hypotheses stated in Proposition 3.11, equation (3.19), are satisfied. 


5. Appendix A 


In this Appendix we collect some useful results obtained in m- Let us consider the 
initial value problem 


(5.1) 


dtu = A{t)u + /(t), n(s) = uq 0 <s<t<T. 


Definition 5.1. Let X a Banach space. A/amily {A(t)}jg[o,T] of infinitesimal generators 
of Cq semigroups on X is called stable if there are constants M > 1 and oj (called the 
stability constants) such that 

(a;,+oo) C p{A{t)), for t G [0, T] 


and 


: A{tj)) < M(A - w) 


-k 


for X > Lo and for every finite sequence 0 < < t 2 , ■■■■,tk <T, k = 1,2,. 


If, for t G [0,T], A{t) is the infinitesimal generator of a Cq semigroup St{s), s > 0 
satisfying ||5t(s)|| < e‘^®, then the family {A(t)}ig[o^T] is clearly stable with constants 
AI = I and to. Precisely, if the operator A{t) generates a Cq semigroup St{s) for every 
fixed t G [0,T], and we can find an estimate for ||5t(s)|| that is independent of t, then the 
whole family {A(t)}^g[Q is stable in the sense of Definition 

Theorem 5.2. Let {A(t)}^g[o,r] be a stable family of infinitesimal generators with stability 
constants M and oj. Let B{t), 0 <t <T be a bounded linear operators on X. // ||il(t)|| < 
K for all t < T, then {A{t) is a stable family of infinitesimal generators with 

stability constants M and oj + MK. 


5.1 
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Now we prove the existence of the so called evolution system U{t, s) for the initial value 
problem (5.1), that is a generalization of the semigroup generated by a linear operator A, 
when such operator depends on time. To this aim, let us state the following result (for 
more details, see m Theorem 2.3, Theorem 3.1, Theorem 4.2]). 


Theorem 5.3. Let {^(t)} 4 e[o,T] a stable family of infinitesimal generators of Cq semi¬ 
groups on X. If D{A(t)) = D, that is the domain of A{t) is independent on t, and for 
uq G D, A{t)uo is continuously differentiable in X, then there exists a unique evolution 
system U{t, s), 0 < s < t < T, satisfying 


\\U{t,s)\\ < for 0 < s < f < T. 


Morevoer, if f ^ C'([s,r],X), then, for every uq G X, the initial value problem (5.1) has 
a unique solution given by 


u{t) = U{t,s)uo-\- U{t,r)f{r) dr. 


for all 0 < s <t <T. 
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